Josephson Current through Charge Density Waves 
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The effect of the colfective charge density wave (CDW) 
motion on the Josephson current in a superconductor/charge 
density wave/superconductor junction is studied theoreti- 
cally. By deriving the kinetic equations for the coupled 
superconductor-CDW system, it is shown that below the criti- 
cal current the CDW does not move. Biased above this value, 
the Josephson current oscillates as a function of the veloc- 
ity of the sliding CDW and the collective mode acts as a 
non-linear shunting resistor parallel to the Josephson chan- 
nel. The effects of impurity pinning are accounted for at a 
phenomenological level 

PACS numbers: 74.50, 72.15.N 



The Josephson effect [Q is known to exist in supercon- 
ductor hybrid structures, where two superconductors are 
separated by insulating barriers or normal metals 0] . Re- 
cently, the theory of the Josephson effect through systems 
which support a non-Fermi (Tomonaga-Luttinger) liquid 
ground state has received much attention due to rapid de- 
velopments in fabrication technology 0. In this Letter 
we investigate the Josephson current through a different 
yet related system, namely a strongly anisotropic metal 
with a charge density wave (CDW) instability. Through- 
out this paper our attention will be restricted to incom- 
mensurate CDW's. Despite the insulating quasi-particle 
spectrum, incommensurate CDW's allow for a unique col- 
lective mode of transport, which distinguishes them from 
ordinary insulators Q . The sliding motion of the CDW 
carries electrical current and is dissipationless in the ab- 
sence of damping sources. It is therefore interesting to 
investigate whether the Josephson current in a super- 
conductor/charge density wave/superconductor (S/C/S) 
junction will be affected by the sliding motion. Recent 
progress in controlled deposition of thin films of CDW's 
may lead to the fabrication of such mesoscopic-scale het- 
erostructures in the near future pj. 

Our model is a current biased S/C/S junction, which 
consists of parallel one-dimensional CDW chains of 
length L, sandwiched between two large superconductors 
with phases ipR and tp L (Fig. 1). The CDW is character- 
ized by a complex order parameter |A c |e lx , where |A C | 
is proportional to the amplitude of the density modula- 
tion and the phase x determines its position. The CDW 
can slide along the chains with velocity proportional to 
X = dtX- 

Using the Keldysh formalism for superconductors ra] 



and CDW's § we have formulated a consistent frame- 
work for the quasi-classical dynamics of the coupled 
superconductor-CDW system. It is shown that the CDW 
is immobile under a certain critical current. Biased 
above this value, the CDW starts to slide and the col- 
lective transport mode causes a shunting resistance par- 
allel to the Josephson channel. We show that, in the 
quasi-stationary approximation, the Josephson current 
has an oscillatory behavior as a function of the phase 
= ( Pr~ ( Pl+xI j / v f, where vf is the Fermi velocity. Ap- 
parently CDW motion induces a dynamical phase which 
is added to the conventional superconductor phase differ- 
ence. Finally, we discuss the effects of impurity pinning 
at a phenomenological level. 

The dynamics of superconductor and CDW systems 
can be described simultaneously by the semiclassical 
Green functions g l a/3 (x;t,t') where i = {R, A, K} and 
a, (3 = {1,2,3,4}. The retarded g R and advanced func- 
tions g A determine the excitation spectrum and the 
Keldysh function g K describes the kinetics of the sys- 
tem. The subscripts 1,3 refer to right and left moving 
electrons with spin up and 2, 4 to left and right moving 
holes with spin down. Throughout this paper "hat" de- 
notes (2x2) matrices and bold-face (4x4) matrices. The 
Green functions satisfy the equation of motion 



ihv F d g g* + H°g 
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Here $ is the quasi-particle potential, a k with k = 
{1,2,3} are the three Pauli matrices. The dot opera- 
tion o denotes internal time integrations as well as ma- 
trix multiplication and the brackets []_ denote commu- 
tation. The self-energy term for impurity scattering is 
neglected throughout this paper. The matrix A s is given 
by A ajU = A s ,22 = 0, A s ,i2 = -A* 21 = |A S | exp(iip), 
and A c = 1|A C | exp(ix). 

It is convenient to gauge away both phases tp and x by 
applying the unitary transformation 
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g* = Ut(x J *)og i oU(a:,t') 



(3) 



where U = exp(|S3X + |<7"3^). Disregarding local vari- 
ations of Xi we look for a stationary state solution of 
the form g(x,t — t), which can be treated by the Fourier 
transformation 



■ie(t—t')/H 



(4) 



The stationary-state equation of motion for the Fourier 
transformed function is 



iv F d x i 



ecr 3 T, 3 - $S 3 - ev F Acr 3 + i\A\, g J 



0. 



(5) 



with $ = $ + \Tiv F d x \ + ^hip, A = Hx/2evF + hd x (p/2e 
and A = |A s |cr2S3 — |A c |S2- From the structure of this 
equation the well-known duality between the supercon- 
ductor and CDW phases Vpd x x ^ <P an d X v Fd x f 
is observed. The gradient of the superconductor phase 
and the time derivative of the CDW phase correspond to 
an electrical current, whereas the gradient of the CDW 
phase and the time derivative of the superconductor 
phase correspond to an electronic potential. 

For the inhomogeneous S/C/S system Eq. (|) has to be 
supplemented by boundary conditions, which adequately 
describe the two superconductor interfaces. Here we re- 
strict ourselves to the ideal case where no defects or po- 
tential barriers are present at the interaces. In order 
to derive the boundary conditions it will be convenient 
to decompose the (4x4) Green functions into four (2x2) 
blocks as follows 



(6) 



It can be shown (Hj] that the diagonal blocks are normal- 
ized as usual 




9° 9 = 9° 9 = lS(t-t'), 



(7) 



and the non-diagonal blocks must satisfy the following 
relations in the superconducting leads on the left (L) and 
right (R) 



9l° }l = -h, 9l ° II = f 



L • 



9r °Ir = Ir, 9r ° fi 



-fi 



(8a) 



(8b) 



The retarded and advanced Green functions in the su- 
perconductors are determined from the stationary-state 
equation of motion Eq. (||) with $ = |A C | = 0. As- 
suming A s to be constant in the superconductors, it is 
convenient to apply the Bogoliubov transformation: 



d+ 




u± —v± 
v± —u± 



(9) 



where u± and v± are the gauge transformed BCS coher- 
ence factors given by 



u± 



A± = sign( e± W4 - |A,P 9(|e ± | - |A a |) 




46(|A S |-| £± |). 



(10) 



Here e± = e^fTix/2, and O is the Heavyside step- function. 
After substitution of Eq. (||), it follows together with 
conditions (0) that the retarded Green functions Q and 
Q at the interfaces are 



Ql = 



Qr = 



1 
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(11a) 
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where the constants Cr^ and Cr l denote Andreev scat- 
tering amplitudes for the scattering of an quasi-electron 
into a quasi-hole and vice- versa. The T and T blocks are 
given by 



Tl 1 Dm D L2 I ' Tl 



Tr 



Dri Dr 2 





Tr 



Dli D L2 



Dri Dr 2 



(11c) 
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where the diagonal D-constants are 'normal' backscat- 
tering amplitudes without spin- flip ('normal' means non- 
Andreev), and the non-diagonal _D-constants are scat- 
tering amplitudes with spin-flip, which will turn out to 
be zero. The boundary conditions ( [ll"| ) state that both 
quasi-electrons and quasi-holes which move away from 
the CDW region into the ideal superconducting leads will 
never be reflected into quasi-particles moving in the op- 
posite direction The boundary conditions for the 
advanced Green functions are obtained by the relation 

S A = -<T3( g fl )t<7 3 . 

In order to determine the yet unknown constants Cs 
& D's, we have to solve Eq. (||) in the CDW region 
and match the boundary conditions. Neglecting the spa- 
tial variations of |A C | near the contacts, the solution 
to Eq. (|) in the CDW region in terms of Bogoliubov 
transformed functions is 



Q{x) = Mr l Q{0)M., 
iRx , 



M = $ exp( 



TiV) 



(12) 



2 



where H = e<r 3 £ 3 — ^x<r 3 — i|A c |S 2 . Evaluating Q{x ) 
at x — L and inserting the boundary conditions (11a)- 
( |lloj ) we obtain the relation which connects the Green 
functions at the left and right interface 



in terms of the transfer matrix C = M.{L): 



C = tf 



exp(-^ 



'<x 3 ) 



cos(^)l 
nvF 



i—sm{- — ) 
A c nvF 



(13) 



H 



«|A C |S 2 



+;^A 



A c p 6(|e| - |A C |) 
^F^e(|A c |-|e|), 



(14) 



with 9 = ipji — ifL + xL/vf- Evidently, the extra phase 
factor arises from the line integral of the vector poten- 
tial A along the junction ^ dxA — xL/vp- Since for 
A c | =/= 0, the solutions are not invariant under the gauge 
transformation, it is expected from Eq. (jl4|) that in ad- 
dition to ipn — (fiL the dynamcial phase xLI v f due to 
the motion of the CDW will also enter the expression 
for the Josephson current. Inserting the boundary con- 
ditions, which describe the two superconductors, we can 
solve for the scattering constants C's and D's. Details 
will be given elsewhere |13| . We next substitute the ob- 
tained solutions into the expression for the total current 
/ through the S/C/S junction 



-1 = 



ev F N(0) 



J deTv(T 3 i K 



ev F N(0)h X . (15) 



with — e the electron charge and iV(0) = (ttTivf)^ 1 the 
density of states at the Fermi level for one spin direction. 
The result is 



I = —X 



-Im / de 



smrj 



, : h J cos r\ — cos £ 



In this expression h q = tanh(^^j^) + tanh( 2 , , 
and we have defined the complex angles 



(16) 



fix/2 n 



2 £± + A± 



(17) 



and 



cos C = 1 — 

' \ C L . 
2 I cos - — sm — 

nvp 2 



e . X C L C+ + 
— sm - — cos — 
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(18) 



In the static limit where x = 0, equation ( |16[ ) corresponds 
to the Josephson current through a band-insulator with 
an energy gap 2|A C |. Taking the limit |A C | — > re- 
produces the well-known expressions for an ideal S/N/S 



junction |l0| , |l2| ]. In the range Tlx < I A* I < |A C | charac- 
teristic oscillations are expected in the Josephson current 
as a function of the CDW velocity. 

Equation ( |l6| ) expresses the total current through the 
S/C/S junction as a function of the superconductor 
phases 4>r,l and the sliding velocity x- In order to obtain 
a closed set of equations, however, it is necessary to derive 
an additional relationship between x and the supercon- 
ductor phase difference. The additional closing relation 
can be derived microscopically from the self-consistency 
relation for the phase of the CDW order parameter 



rfeTrS 3 g A ' = 



(19) 



and the Keldysh functions in the reservoirs. The self- 
consistent solution for the clean junction requires an elec- 
trochemical potential difference 5/i = (J,r — fJ>L in the 
sliding state Sfi — fix, as in Ref. ||. This implies that 
below the critical Josephson current the CDW will not 
move. Biased above the critical current the CDW slides, 
but since the collective CDW motion is dissipative due 
to the contact reservoirs, a small potential difference is 
induced on the superconducting leads and the supercon- 
ductor phase-difference will evolve slowly in time. As a 
consequence we obtain the relation 



Sfi = hip/2 = hx- 



(20) 



The sliding CDW mode thus acts as a shunting resis- 
tor parallel to the Josephson channel. The above quasi- 
stationary approximation, where the dynamics of the su- 
perconductor phase is neglected on the Josephson term, 
is therefore justified a posteriori pT[ . In the short junc- 
tion limit L <C 7u^f/|A s |, the dynamical phase x^/vf 
can be neglected and if |A C | ^> |A S | the Josephson cur- 
rent can be written as I c sin ip, where I c is the critical 
current. The time-averaged potential becomes < V >— 
Rq ^Jl 2 — I 2 and it oscillates with the Josephson fre- 
quency L) j = 2e < V > /h, where Rq = h/2e 2 is the 
(spin included) quantum resistance of a one-dimensional 
ballistic quantum wire. Note that far above the crit- 
ical current I I c the Josephson frequency becomes 
two times the fundamental frequency of the narrow band 
noise wj — 2ljnbn — 2irl/e. 

The assumption of an ideal S/C/S junction is rather 
restrictive in realistic systems where interface effects and 
impurities tend to pin the CDW. Although a fully mi- 
croscopic treatment is beyond present work, one can ac- 
count for the pinning effects by replacing Eq. (|2(]) by 
5[i = Tup/2 = fix + /J-t sin x- This model corresponds 
to the single-particle model of CDW's 14 1, where /it 
represents the threshold value required to overcome the 
pinning potential. 

We conclude by summarizing our results. We 
have formulated the kinetic equations for a coupled 
superconductor-CDW system and calculated the cur- 
rent through an ideal superconductor/charge density 



3 



wave/superconductor junction. In the dc limit the CDW 
does not move. Biased above the critical current the 
CDW slides, and the Josephson current oscillates as a 
function of the CDW velocity. The collective mode acts 
as a non- linear shunting resistance parallel to the Joseph- 
son channel. 
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FIG. 1. Schematic figure of a charge density wave with 
length L, sandwiched between two superconductors with 
phases <f>n and <f>L- The junction is biased with a current 
source. 
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